Hofstadter butterfly as Quantum phase diagram 
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The Hofstadter butterfly is viewed as a quantum phase diagram with infinitely many phases, 
labeled by their (integer) Hall conductance, and a fractal structure. We describe various properties 
of this phase diagram: We establish Gibbs phase rules; count the number of components of each 
phase, and characterize the set of multiple phase coexistence. 
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Introduction. — Azbel Q recognized that the spectral 
properties of two-dimensional, periodic, quantum sys- 
tems have sensitive dependence on the magnetic flux 
through a unit cell. A simple model conceived by Peierls 
and put to the eponymous Harper as a thesis problem, 
gained popularity with D. Hofstadter Ph.D. thesis |Q, 
where a wonderful diagram, reminiscent of a fractal but- 
terfly, provided a source of inspiration and a tool for spec- 
tral analysis 

The Hofstadter butterfly can also be viewed as the 
quantum (zero temperature) phase diagram for the in- 
teger quantum Hall effect. It is a fractal phase diagram 
with infinitely many phases jLiyiSj ]. The diagram leads 
to certain natural questions: Count the number of com- 
ponents of a given phase; Classify which phases coexist 
and where. It also leads to the general question: What 
form does the Gibbs phase rule take for quantum 

phase transitions. 

Fractal phase diagrams and/or infinitely many phases 
appear in dynamical systems 
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In classical lat- 
tice systems fractal phase diagrams |l4| , p^| , po| , p^| are com- 
monly viewed as a pathology due to either long range in- 
teractions, or, as is the case for spin glasses, loss of trans- 
lation invariance. The Hofstadter model, when viewed as 
a statistical mechanical model, is both short range and 
translation invariant in a natural way. But, it is quan- 
tum and the translation group is non-commutative. It 
suggests that fractal phase diagram may be more com- 
mon in quantum phase transitions than in classical phase 
transitions. 

The Hofstadter Model. — The model conceived by 
Peierls has two versions. For the sake of concreteness 
we shall focus here on the tight binding version. On the 
lattice Z 2 , define magnetic shifts 



(Uif))(n, to) = ip(n — 1, m) 
(V($)V0 (n, to) = e 27rin *^(n,TO - 1) 
n, to £ Z 



(1) 



2irQ is the magnetic flux through a unit cell. The Hofs- 
tadter model is 

H ■(*, a, b) = a(U + U*) + &(V($) + V* ($)), (2) 

where a, b > are "hopping" amplitudes, a — b is called 



the self-dual case |1C|] and we shall focus on that case 
below. We set = if ($, 1, 1). 

U and V, and therefore also H($), commute with the 
(dual) magnetic translations U and V, 



(U%p){n, m) = ip{n, m — 1) 



(Vijj)(n, to) 



ip{n — 1, to) 



(3) 



This makes H{<&) translation invariant in a natural 
way. The group of magnetic translations pi] ] is non- 
commutative: 



U*V*UV = UVWV* = e~ 27ri *. 



(4) 



The one particle representation of the Hofstadter model, 
Eq. (0), is natural for spectral studies. In the context 
of statistical mechanics the second quantized representa- 
tion of the model is also instructive because it makes it 
clear that the model has short range, in fact, only on site 
and nearest neighbors, interactions. The fractal features 
of the phase diagram are, therefore, not a consequence of 
long range forces, as in some classical statistical mechan- 
ics models. The second quantized form is: 



W(*,m) = E ei7(nm; "' m ' )a » m <W 
where 

f 1 



[i > a! 



(5) 



n — n' = ±1, m = to'; 
to — m' = ±1, n = n'; (6) 
otherwise. 



fi is the chemical potential and , a are the usual 
Fermionic operators. 

Let us recall few elementary features of the spectrum: 



S(H($)) = -S(H(&)) = -S(H(1 - $)). 



(7) 



The first is a consequence of Z 2 being bipartite, and the 
second is a consequence of time reversal. Together, they 
imply a four-fold symmetry, manifest in the Hofstadter 
butterfly. 

The electronic density, p(<I>, fx), (=integrated density of 
state), is 



(8) 
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where |0) is Kroneker delta at the origin. < p(<I>, p) < 1, 
is an increasing function of p. 9 is the usual step function. 

The gaps in the spectrum are labeled by an integer, k, 
which is a solution of p5p| . 



where <f>(j) is Euler (totient) function. Recall that (f>(j) 
counts the number of integers, up to j (and including 1), 
that are prime to j: (f)(1) = 1, (f)(2) = 1, (f>(3) = 2 etc.. 
To prove Eq. (12) note first that from Eq. (||) 



$ k — p mod 1 . 



(9) 



P = P 



(13) 



k is the Hall conductance. We picked the letter k be- 
cause it is naturally associated with an integer, and it is 
also the first letter in von Klitzing name. By Eq.(^) and 
Eq. © : 

k(p,$) = -k(p,l-<S>)=-k(-p,<f>). (10) 

which implies a four fold (anti) symmetry of the butterfly. 

We shall assume that the Ten Martini conjecture p6| ] 
holds. Namely, that for all irrational <i>'s, all the gaps are 
open, so Eq. (||) has p in an open gap for all k e Z. 

Fig. [j] shows the Hofstadter butterfly, color coded ac- 
cording to the Hall conductance. The color coding counts 
the Hall conductance modulo 16. Zero Hall conductance 
and the spectrum are left blank. The gross features of 
the diagram are associated with small integers where the 
color coding is faithful. 

The colored picture emphasizes the gaps while the 
standard Hofstadter butterfly emphasis the spectrum. 
The colored figure is prettier and displays the regular 
aspects of the diagram: Gaps are better behaved than 
spectra. The colored diagram is also more faithful to cer- 
tain spectral characteristics. For example, the spectrum 
is a small set, (in fact, one of zero Lebesgue measure), 
something that is manifest in the colored diagram, but 
is less obvious from the usual Hofstadter butterfly which 
plots the spectrum. 

We also broke with tradition in that the colored Hofs- 
tadter butterfly is rotated by 90°: In Fig. | the horizontal 
axis is $ and the vertical axis is the energy, or S(H (&)*). 
The reason we chose to do so is that this way emphasizes 
the fact that phase boundaries are functions (of $). 

We denote by P(k) the fc-th phase. Formally, 

P(k) = {$, p\ $fc = p($, p) mod 1, p i S(H(<f>)) } (11) 

P(k) is an open set in the (<&, p) plane, with a finite num- 
ber of components. For example, P(l) is two of the four 
big wings of the butterfly. We call P(k) a pure phase and 
denote its number of components |P(fc)|. The closure of 
the pure phase is denoted P(k) and the phase boundary 
is dP(k). We call dP(k) the total boundary. 

Counting components. — The A:-th pure phase is made 
of several components. The k = phase (blank) has two 
components. For k the number of components is 

2|fe| 2 
|PW|^,t(j) = 12- + 0(Ho g i) (12) 



Hence, a given color would appear at most once on any 
vertical line of fixed $. 




FIG. 1. Hofstadter colored butterfly 



When $ = | with gcd(p, q) = 1 (gcd(p, q) is the great- 
est common divisor of p and q) , the spectrum of the Hof- 
stadter model has q finite bands and q — 1 gaps which 
are all open intervals (except when q is even the central 
gap is closed) . We then number the gaps by their natural 
order 1, . . . , q — 1. The semi-infinite interval below the 
spectrum is, formally, the 0-th gap, and the semi-infinite 
interval above the spectrum as the g-th gap. Eq. (^), for 
the j-th gap, can be written as 

pk = j mod q. (14) 

Given p, q and j the equation has a unique solution for 
each open gap such that |fe| < S. In particular, once q 
has been fixed, the Hall conductance takes all (non-zero) 
integer values, from — [|J to [_§ J , and each value appears 
once. For even q, the central gap is closed and formally 
can be assigned a value of ±|. 
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The number of components of the A;-th phase is the 
same as the number of flux values, which accommo- 
date the wing-tips, minus one. The fc-th wings-tips are 
located at those values of <!> where the fc-th color is ab- 
sent (the fc-th color is present in any small neighborhood 
of those values of <&, because this neighborhood contains 
fractions with arbitrary large g's). Given k ^ 0, it must 
appear once on a horizontal interval with = ^ pro- 
vided | > |fc|. The tips of a wing with a given color 
must, therefore, be located at those values of <& which do 



not admit A; as solution of Eq. (14) for any p. In other 
words, the wing-tips lie at those values of <& where q is 
too small to accommodate \k\. This is the finite set, a 
Farey sequence, 

F m = | 0<p<q,gcd(p,q) = l,q<2\k\} (15) 

2|fc| . . 

= Ui °<P< 1>9cd(p,q) = l\\J{0} 

q=l ' 



Let \F\ be the number of elements in F. Then from (15) 



\P(k)\ = \Fr 



2\k\\ 



2\k\ 



(16) 



9=1 



essentially from the definition of the Euler function. 

The asymptotic expansion for the sum in Eq. (|l2|) is 
taken from |2^] . 

Pure phases and Phase boundaries. — In thermody- 
namics and statistical mechanics, Gibbs phase rule is a 
statement about the structure of pure phases and their 
boundaries. A weak form of the Gibbs phase rule says 
that pure phases are a set of full measure; two phases co- 
exist, generically, on a set of Hausdorff co-dimension one 
etc. [Q. This is the form that one gets if one considers 
general convex functions for thermodynamic potentials. 
The number of coexisting phases is related to the dimen- 
sion of tangent planes, and the Gibbs phase rule is a con- 
sequence of theorems about convex functions. There is a 
stronger form of the rule Jig] ] which posits, in addition, 
that the sets are (locally) manifolds. This form is a con- 
sequence of additional regularity of the thermodynamic 
potentials. 

Gibbs phase rule is a consequence of the convexity of 
thermodynamic potentials and so is ultimately based on 
the second law of thermodynamics. It has nothing to 
say about the zero temperature phase diagram of quan- 
tum phase transitions in general, |]l5]| , and the Hofstadtcr 
model in particular (because the entropy vanishes iden- 
tically). A question that arises is then what form might 
Gibbs phase rule take for quantum phase transition. The 
Hofstadtcr model restricts what could and what could 
not be true in general. As we shall see, the phase dia- 
gram of the Hofstadter model turns out to satisfy only a 
weak form of the Gibbs phase rules. 



Fig. [j] suggests that the set of unique phase is a set 
of full measure and that the phase boundaries, though 
fractal, are not too wild. More precisely, we have: 

Gibbs-like phase rule: The phase diagram of the 
self- dual Hofstadter model is such that pure phases, 
labeled by their Hall conductances, are full measure; 
phase boundaries are not manifolds — they are nowhere 
differentiable- but they are almost so in the sense that 
their Hausdorff co- dimension is integral, in fact: 



dim H (dP(k)) = 1. 



(17) 



Since the number of phases is countable the total phase 
boundary U m {dP{mf) , is a set of Hausdorff dimension 
one as well. Finally, infinitely many phases coexist on a 
countable set, and therefore a set of Hausdorff dimension 
zero. 

The first part of the Gibbs-like phase rule is an easy 
consequence of a result of Last ||, which states that 
\S(H(§)) J = for a set of <I> of full measure. That phase 
boundaries are nowhere differentiable follows from results 
of Wilkinson j9), Rammal, and Helffer and Sjostrand j|] 
who showed that the phase boundaries dP(k) possess dis- 
tinct left and right tangent at every rational $. That the 
Hausdorff dimension of the boundary is one follows from 
results of Bellissard jl6| who showed that away from the 
wings tips, dP(k) can be represented by a functions of 
$ that are uniformly Lipshitz. By standard results, p5| ] 
it then follows that the Hausdorff dimension is one. The 
set of infinite phase coexistence is analyzed below. 

Coexistence. — In [|l8| the term lakes of Wada was used 
to describe dynamical systems with the property that any 
point on the boundary of the one basin of attraction is 
also on the boundary of all other basins. We shall say 
that a system is almost Wada of order m if every circle 
that contains two pure phases contains m pure phases. 

The Hofstadter butterfly is almost Wada of infinite or- 
der. This is seen from the figure, and can also be shown 
to follow from Eq. (||). 

We say that the two pure phases, P(m), P(n), coexist 
on 



C(m,n) = dP{m) n dP(n) 



(18) 



No two phase coexist for any irrational flux. This is easily 
seen from Eq. ([)]): For irrational $ the electron density p 
takes a dense set of values in the gaps. Therefore, any two 
phases, P(m) and P(n) are separated by infinitely many 
other phases. It follows that the set of phase coexistence 
is a countable set, and so of zero Hausdorff dimension. 

The following result gives a complete characterization 
of phase coexistence: 

Proposition: Consider a point x <E dP{k) with 
$(x) = H with gcd(p,q) = 1. Then x G dP(k + Iq) 
for all £ £ Z. Moreover x £ dP(k') if k' ^ k + iq for 
each £ GZ 
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Proof: Since gcd(p, q) = 1 the equation 
pa — qb = 1 



(19) 



has a solution with integer a and b (where a is non-unique 
mod q). Let ^ = with neZ. Then 



PnQ- QnP = 1. 



(20) 



From Eq. ( [H] ) it follows that each band at ^2. carries Hall 
conductance g (mod g n ). 

Now consider a point a; on the right boundary of P(k). 
We shall first show that x is also on the left boundary 
of P(k + q). Since P(k + q) has a finite set of wing- 
tips, Eq. (|15|), when n is large enough, the gap with label 
(=Hall conductance) k + q at flux |a must be open, and 
must remain open for all large n. By a bound of Last 
and Wilkinson [E3 for the total width of the spectrum 



at — , each band is small and hence 



24 

dist(P(k),P(k + q)) < —. 



(21) 



Taking n — > 00 we see that a; € 0P(fc + q) as claimed. 

By considering the next band we shall now show that x 
also lies on the boundary of P(k + 2q). Now, P(k + 2q) at 
Y~ is separated from P(k) by two bands and a gap. The 
bands are small by the Last Wilkinson bound. The gap 
is also small by the Holder continuity of the spectrum: 



\gap\ < 18 



Pn 

q n 



p 
<1 
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and from this 



dist{P(k),P(k + 2q)) < 



24 



18 



(23) 



Taking the limit n — > 00 yield the result. The argument 
can be repeated for any P(k + tq) with I finite and pos- 
itive. Negative values P(k — lq) are obtained by letting 
n — > — 00 in the argument above. 

For the left boundary point of P(k), n — > 00 will give 
P(k — lq) and n — ► —00 will give P(k + lq). This for- 
mula applies also to the phase k = 0, with it's right 
boundary being the leftmost point of the spectrum and 
vice- versa. The Hall conductance for the middle gap with 
even q (which is closed) is formally ±|, so it is common 
to phases P(±§ +l'q) which is the same as P(| + lq) for 

e, e el 

The second part of the proposition follows from the 
equality 



U Ike*- 



lq) = |J P(k) 



fcez 



and the fact the for fixed <i>, each Hall conductance k can 
appear only once. 
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